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Проведено численное моделирование движения вязкой несжимаемой неизотермической жидкости в открытой полости прямоугольного 
типа в условиях вынужденной конвекции и сопряжённого теплообмена. Получена гидродинамическая картина течения вязкой жидкости 
в открытой полости в сопряжённой и несопряжённой постановках. Получены температурные профили для двух фаз - твёрдой и жидкой. 
Изучено влияние параметров модели на характер движения. Показано влияние параметров модели на характер распределения температуры 
в обеих фазах. 


Устойчивый интерес к исследованиям конвективных 
течений в полостях различных форм и типов наблюда- 
ется на протяжении последних сорока лет. Интерес этот 
объясняется широким прикладным значением пробле- 
мы: полости в качестве теплопередающих, теплоизо- 
лирующих и технологических элементов встречаются 
в энергетических и технологических установках различ- 
ного предназначения, радиоэлектронных устройствах и 
теплообменной аппаратуре [1,2]. 

Исследование лобового взаимодействия струи вяз- 
кой несжимаемой неизотермической жидкости (ВННЖ) 
с ограниченным объёмом различной формы имеет 
важное научно-практическое значение в связи с тем, 
что подобные течения широко распространены в тех- 
нологических процессах различного уровня сложности 
таких отраслей промышленности как энергетическая, 
нефтехимическая, атомная, металлургическая, аэрокос- 
мическая и многих других [2-5]. 

Моделирование теплообмена при движении вязкой 
жидкости в полости прямоугольного типа сопряжено 
с решением достаточно сложных задач вынужденной 
конвекции несжимаемой жидкости. Поскольку создание 
надёжных аналитических методов расчёта параметров 
течения ВННЖ в ограниченных объёмах различного 
типа исключено из-за сложности таких течений, то 
возникает необходимость численного моделирования. 

В данной работе рассматривается нестационарное 
взаимодействие ламинарной струи ВННЖ с открытой 
полостью прямоугольного типа (рис. 1). Цель данной 
работы - исследовать гидродинамику и сопряжённый 
теплообмен при движении ВННЖ в открытой полости 
прямоугольного типа в условиях вынужденной кон- 
векции. 

Изучение описанного процесса проводилось с 


использованием математической модели на основе 
системы уравнений Навье-Стокса в переменных вихрь 
- гЪѵгплшя тпк-я при умеренных числах Рейнольдса 
1 00 < Ке < 800 ) уравнения энергии, а также уравнения 
теплопроводности для материала прямоугольной по- 
лости с соответствующими начальными и граничными 
условиями: 
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Численное решение гидродинамической задачи 
осуществлялось в области - 3, ограниченной участком 
затекания - 1, линией симметрии -4, боковой стенкой 
- 6 и дном полости - 5, а также участком выхода - 2(3] 
прямоугольной ПОЛОС'”" 1 ' 

На дне полости ѵУ ~ з,1Д<х-. ) высхавляехся 
условие непротекания, прилипания, а также граничное 
условие четвёртого рода для уравнения энергии: 
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Условия, аналогичные описанным выше - условие 
непротекания, прилипания, а также граничное условие 
четвёртого рода для уравнен™ ™ тган ^тг'г авляются 
на боковой стенке полости ( х = ІГ8<у<Н) : (З) 

ѵ|/ = 0; от = 2-хр(х + Ах, у)/(Ах) 2 ; 
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Условия неразрывности тепловых потоков и 
непротекания заданы на оси симметрии струи ( 
х=Ь, 8<у<Н у 
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При затекании жидкости в полость, как уже упо- 
миналось, выделяются два участка - участок входа в 
полость - 1 и участок выхода из полости - 2 (рис. 1). 
При исследовании данного процесса необходимо вы- 
полнение интегрального соотношения, определяющего 
расход жидкости: 

( 8 ) 
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На участке выхода из полости прямоугольного типа 
для составляющих скорости использовалось условие 
“сноса”: 
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для температуры мягкое граничное условие - вторая 
производная температуры по координате у равна О 
[6-8]. Условия теплоизоляции задаются на внешних 
границах полости 
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ем движения жидкости в полости ( 1 < 

- координата фиксированной крайней точки 
участка, лежащей на боковой стенке полости; ѵ — 1 х - ^ ' 

- поперечная составляющая скорости движения жид- 

кости ] у / х ѵ \нии от входного участка до дна по- 
лости; + " - поперечная составляющая скорости 

движения жидкости в направлении от дна полости до 
участка выхода. 
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Рис. 1 . Общая схема течения в прямоугольной полости и геометрия 
расчётной области: 1) участок входа в полость; 2) участок 
выхода из полости; 3) гидродинамическая область; 4) ось 
симметрии; 5) дно полости; 6) боковая стенка полости; 7) 
внешняя стенка полости 


Методом конечных разностей решены система 
уравнений Навье-Стокса в переменных функция тока 
- вихрь, уравнение энергии и уравнение теплопро- 
водности для (1-12). Разностные аналоги уравнений 
переноса и теплопроводности решены методом про- 
гонки [7-9]. Уравнение Пуассона на каждом временном 
слое решалось методом последовательной верхней 
релаксации. Использовалась разностная схема вто^^] 
порядка точности. Расчёты выполнялись на равномер- 
ной разностной сетке. 

При постановке задачи для системы уравнений 
Навье-Стокса в переменных вихрь - функция тока 
главной особенностью являются граничные условия, 
которые в случае твёрдой неподвижной поверхности 
имеют следующий вид [7, 8]: 

,=0, *-0. 

Здесь п - нормаль к твёрдой поверхности. 

При численном решении разностного уравнения 
для вихря возникает проблема определения граничных 
условий для него, поскольку представленные граничные 
условия (13), относящиеся к системе уравнений Навье- 
Стокса, заданы только для функции тока и не заданы для 
вихря. В данной работе использовался способ, заключа- 
ющийся в том, что функция тока вблизи границы пред- 
ставляется в виде ряда Тейлора. Выражение для вихря 
на границе получаем после проведения необходимых 
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операций с разложением функции тока [6]. 

Для проверки аппроксимации и сходимости числен- 
ного решения в качестве тестовой была взята задача о 
плоском течении в прямоугольной каверне с верхней 
стенкой, движущейся в своей плоскости с постоянной 
скоростью [6 01 Сравнение резул т ' го тов по профилям 
продольной и ' - * и поперечной ѵ ' У' составляющих 
скорости с данными других авторов [6-8], показало 
хорошее, в пределах ±7 %, согласование, а также по 
значениям Лшнкчши тпк-я при различных числах Рей- 
нольдса -ТОО — іѵ С — 1 000 „ попптііітглл| ѵлтігтл^па 

узлов разностной сетки 20 < к < 60, 17<]<60, 
где к, ] - количество узлов сетки по координатам х, у 
соответственно. 

В качестве второй тестовой была решена задача о 
сдвиговом течении г тшпк-уляцией при малых числах 
Рейнольдса Ю<К.е<50 щд] Решалось бигармони- 
ческое уравнение для функции тока, так как известные 
подстановки исключают функцию вихря и непосредс- 
твенно определяются поля неизвестной функции тока 
[10-12]. Сравнение результатов по профилям функции 
тока в различных сечениях с данными, полученными в 
монографии [10], показало хорошее, в пределах ±5 %, 
согласование. 

В ходе проведения научно-исследовательской ра- 
боты рассматривались жидкости разных типов (вода, 
расплавленный свинец, жидкая сталь) с широким 
диапазоном изменения динамического параметра Ке и 
параметров модели. В данной статье представлены ре- 
зультаты, полученные при исследовании гидродинамики 
и сопряжённого теплообмена при движении жидкой ста- 
ли в полости, которые имеют большое научно-практи- 
ческое значение как из-за возможности возникновения 
аварийных ситуаций, так и в связи с внедрением новых 
современных технологий и усовершенствованием уже 
существующих в таких отраслях промышленности как 
металлургическая, энергетическая и многих других. 
Температура натекающей жидкости (жидкой стали) 
составляет 1 500 °С. На рис. 2-8 приведены типичные ре- 
зультаты численных исследований для жидкой стали. 

Как следует из анализа установившегося поля 
течения, которое изучалось для различных вариантов 
геометрических характеристик полости (в частности, 
для Ь/Н=1/2, 2/3, 1) при достаточно 
зоне изменения чисел Рейнольдса - ке _ эіііі, 
жидкость доходит до дна выемки, разворачивается и 
вытекает на всём участке 2 (рис. 1 ). Исходя из этого при 
исследовании процесса движения вязкой несжимаемой 
неизотермической жидкости в полости можно выделить 
два этапа. К первому этапу можно отнести прохождение 
жидкости от входного участка до дна полости с учётом 
взаимодействия с ним. Течение при взаимодействии 
струи с дном выемки сопровождается торможением 
жидкости и возникновением области с повышенным 
давлением, что приводит к растеканию жидкости вдоль 
дна полости. Второй этап движения жидкости проходит 
от дна полости до выходного участка, образуя область 
возвратного течения. На данном этапе продолжается 
торможение жидкости, в результате чего возникает так- 
же область с повышенным давлением. Области прямого 


и возвратного течения, соответствующие описанным 
этапам движения жидкости в полости, хорошо видны 
на рис. 4. 

Распределение функции тока х ' - для момента 
времени, соответствующего установившемуся полю 
течения (данный момент времени будет использоваться 
в дальнейшем и для простоты он будет обозначаться в 
тексте как фиксированный момент времени), при числе 
Ке=200 и геометрическом отношении сторон полости 
Ь/Н=1/2 представлено на рис. 2. Максимум функции 
тока соответствует зоне 0,86<у<0,95, в которой наблю- 
дается наиболее интенсивное формирование вихревых 
структур. Необходимо заметить, что проведённый 
анализ позволил выявить достаточно значительное вли- 
яние геометрических характеристик на формирование 
поля функции тока. Изменение входного параметра а 
в сторону его уменьшения позволяет сделать вывод о 
незначительном изменении характера распределения 
изолиний функции тока. 

На рис. 3 показано распределение продольной 
составляющей скорости и(х,у) в фиксированный мо- 
мент времени при числе Ке=200. Максимум скорости 
соответствует зоне 0,40<х<0,44 около дна полости. При 
увеличении числа Рейнольдса характер поведения про- 
дольной составляющей скорости в различных сечениях 
полости качественно сохраняется. 

На рис. 4 показано распределение поперечной со- 
ставляющей скорости ѵ(х,у) в фиксированный момент 
времени при числе Ке=200. Анализируя распределения 
поперечной составляющей скорости, нужно отметить, 
что с увеличением числа Рейнольдса профиль ѵ(х,у) в 
начальных сечениях полости становится более запол- 
ненным и близким к постоянному значению, в то время 
как при числах Рейнольдса Ке=100, 200 профили по- 
перечной составляющей имеют почти параболическое 
распределение в тех же сечениях. Анализ распределения 
поперечной составляющей скорости показывает, что в 
поле течения формируются два максимума, соответс- 
твующих прямому и возвратному течению. На рис. 3, 
4 можно заметить, что на первом этапе по мере про- 
движения струи несжимаемой жидкости к основанию 
выемки поперечная составляющая скорости падает при 
возрастании продольной составляющей, что отчётливо 
наблюдается в области 0,41<у<0,55. На втором этапе 
по мере продвижения струи к выходному участку про- 
дольная составляющая скорости падает - это хорошо 
видно на участке 0,32<х<0,41, 0,41<у<0,55, а поперечная 
составляющая скорости начинает расти - наиболее явно 
это проявляется в области 0,59<у<0,75, 0,30<х<0,40. В 
соответствии с представленным рис. 4 максимальное 
значение скорости вытекания ѵ тах =0,24. 

Рассматривалось влияние на характеристики тече- 
ния длины входного участка а - На рис. 5 приведено 

расп поперечной составляюще” т; ; \ 

а= 0,28а. / N а= 0,16а» 

для ’ “плошные линии) и для 

(пунктирные), а * - обозначение длины проницаемого 
участка полости у=Н, Б<х<Ь. С увеличением длины 
входного участка максимальное значение поперечной 
составляющей скорости достаточно незначительно 
уменьшается на проницаемых участках полости (учас- 
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тки входа и выхода). Характер поведения поперечной 
составляющей скорости в различных сечениях полости 
качественно сохраняется с изменением длины входного 
участка. 

На рис. 6 приведено распределение продольной 
~ 0 2 8 іей скорости для длины участка затекания 

а " а> ' С уменьшением длины входного участка 

максимальное значение продольной скорости снижа- 
ется, и меняется характер распределения функции и 
от х. 

Исследование движения вязкой несжимаемой жид- 
кости в полости прямоугольного типа проводилось в 
условиях сопряжённого теплообмена. Были получены 
температурные зависимости в твёрдой и жидкой фазах 
при различных динамических параметрах и варьирова- 
нии параметра а - 

На рис. 7 представлены изолинии температуры в 



Рис. 2. Распределение функции тока на плоскости х, у в фиксиро- 
ванный момент времени при числе Ке=200 и геометричес- 
ком отношении сторон полости ІѴН=1/2 



плоскости х, у в фиксированный момент времени при числе 


Ке=200 


фиксированный момент времени при Ке=200, Рг=0,979. 
Следует отметить, что при возрастании входного па- 
раметра а и понижении геометрических параметров 
картина поведения температурных полей по высоте 
полости несколько меняется. 

Распределения температуры в твёрдой и жидкой 
фазах в различных сечениях у: 0,36, 0,52, 0,624 в фик- 
і момент времени при Ке=300, Рг=0,979, 
а ’ а * приведены на рис. 8. Как видно из рисунка, 
температурный профиль до границы раздела фаз х=0,25 
убывает достаточно быстро, в то время как в твёрдой 
фазе температура изменяется гораздо медленнее. 

На основании полученных результатов можно 
сделать вывод о том, что характер поведения темпе- 
ратурных профилей в различных сечениях полости 
качественно сохраняется с изменением динамического 
параметра Ке, в то время как входной параметр ока- 
зывает достаточно значительное влияние на характер 
распределения температурных полей и гидродинами- 
ческую картину течения. 

В процессе исследования была также решена не- 
сопряжённая задача. На стенках полости выставлялись 
условия теплоизоляции (отсутствие теплоотдачи в 
стенки полости). 

Сравнение температурных профилей на стенках 
полости, полученных в сопряжённой и несопряжённой 
задачах, показало отличие в пределах ±23 % по значе- 
ниям температуры, что говорит о целесообразности 
сопряжённой постановки. 

Результаты численного анализа позволяют сделать 
вывод о возможности дальнейшего расширения облас- 
ти применения математического аппарата [7, 8] для 
решения задач о конвективных течениях в полостях 
открытого типа в условиях струйного затекания и 
сопряжённого теплообмена. Данная работа является 
логическим продолжением трудов [8, 9], в которых была 
впервые показана возможность применения математи- 
ческого аппарата [7] для решения задач в сопряжённой 
постановке для областей с более сложной геометрией, 



Рис. 4. Распределение поперечной составляющей скорости на 
плоскости х, у в фиксированный момент времени при числе 
Ке=200 
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чем канал или обтекаемое тело [13]. В дальнейшем 
получить устойчивые решения задач для областей со 
значительно более сложной геометрией, чем рассмот- 
ренные в [7-9] и в данной работе, позволит оптимизация 
сеточных параметров в связи с расширением возмож- 
ностей современных ПЭВМ. 

В ходе изучения рассматриваемого процесса специ- 
ально исследовался вопрос об устойчивости численного 
решения рассматриваемой задачи при увеличении числа 
Ке. Согласно проведённому анализу было установлено, 
что устойчивые решения рассматриваемой задачи реа- 
лизуются ° ппртэтшшл тт™ілг, ком диапазоне изменения 
чисел Ке О 60 < Ке < 1 000). 



Рис. 5. " „ іие поперечной составляі 2 , , ти 

а ~ 0 . 28 а. , ч а= 0 , 16 а., 

(сплошные линии) и для (пунк- 

тирные) в сечениях: 1) у=0,675; 2) у=0,567 



Рис. 7. Изолинии температуры в фиксированный момент времѵ..,. 
при К.е=200, Рі=0,979 



Рис. 6. " іие продольной составляющей скорости для 

а ’ а в сечениях: 1) у=0, 324; 

2) у=0,2 16; 3) у=0,135 



Рис. 8. Распределения температуры в твёрдой и жидкой фазах в раз- 

личных сечениях у по координате х в фикси 
мент времени при числах К.е=300, Рг=0, 979, а ' а, ‘ 

1 ) у=0,36; 2) у=0,52; 3) у=0,624; 4) граница раздела жидкой 
и твёрдой фаз 
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